set term jpeg
set output "q5even.jpg"
set xlabel "A(angstrom)"
set ylabel "Xi"
a = 0.3
b = 2.58161/a
f(x) = sin(b * x)
set xrange [-a:a]
set title "Even function with B = ". sprintf("%g",b)
plot f(x) 

set output "q5odd.jpg"
b1 = 1.3000183/a
b2 = 3.8110962/a
f(x) = cos(b1*x)
g(x) = cos(b2*x)
set title "Odd functions"
set xrange [-a:a]
plot f(x) title sprintf("B=%g",b1), \
g(x) title sprintf("B=%g",b2)
